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ABSTRACT 


The structure and dark matter halo core properties of dwarf spheroidal galaxies 
(dSphs) are investigated. A double-isothermal (DIS) model of an isothermal, non self- 
gravitating stellar system, in gravitational equilibrium and embedded in an isothermal 
dark halo core provides an excellent fit to the various observed stellar surface density 
distributions E*(r). Despite its constant veloc ity dispersion, the stellar system can be 
well characterised by King profiles (King 1 9661 ) with a broad distribution of concentra¬ 
tion parameters c = logwith r*^ and r* jC the stellar ti dal and core radius, 
respectively. The DIS model confirms the suggestion of Kormendv &; Freemanl (2014) 
that the core scale length of the stellar system, defined as a* = — (din Xb/dr 2 ) -1 / 2 , is 
sensitive to the central dark matter density pq. In contrast to single-component sys¬ 
tems, however does not trace the tidal radius of the galaxy but the core radius r c 
of its dark matter halo, c is therefore sensitive to the ratio d*/do with d* and do the 
stellar and dark matter velocity dispersion, respectively. Simple empirical relationships 
are derived that allow to calculate the dark halo core parameters po , r c and (To, given 
the observable quantities d*, a* and c. The DIS model is applied to the Milky Way’s 
dSphs. Their halo velocity dispersions lie in a narrow range of 10 km/s < do < 18 
km/s with halo core radii of 280 pc < r c < 1.3 kpc and r c ~ 2a*. All dSphs follow 
closely the same universal scaling relations (por c ) = po x r c = M 0 pc~ 2 and 

do x r c 1 = 0.45^0 27 (km/s) 2 pc -1 that characterise the cores of more massive galaxies 
over a range of 18 magnitudes in blue magnitude Mb- For given (por c ) the core mass is 


a strong function of core radius, M c 


T 

' n 


Inside a fixed radius r u , with r u the logarith¬ 


mic mean of the dSph’s core radii, the total mass M u = 2.17 (po r c) r u is however roughly 
constant. Outliers with smaller masses are expected for dSphs with core radii that are 
much larger or smaller than r u . For the Milky Wa y’s d Sphs we fin d r„ . = 400 ± 100 pc 
and M u = 2.6 ± 1.4 x 10'M 0 , in agreement with Str i gari et al. ( 2008 ). Due to their 
small r c , the core densities of the Galaxy’s dSphs are very higher, with po = 0.03 - 0.3 
M 0 pc~ 3 . The dSphs would have to be on galactic orbits with pericenters smaller than 
a few kpc in order for their stellar systems to be affected by Galactic tides which is 
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very unlikely. dSphs should therefore be tidally undisturbed. Observational evidence 
for tidal effects might then provide a serious challenge for the cold dark matter scenario. 

Subject headings: dark matter - galaxies: dwarf - galaxies:formation - galaxies: kine¬ 
matics - galaxies: structure 


1. Introduction 


The nature of dark matter is still a mystery. Standard cosm ology wo rks with th e as sumption 
of a massive, weakly and g ravitationally interacting particle (e.g. IWhitelIl982l : ISteigman &; Turner 
1985; Jungman et ah 1996j). This cold dark matter (CDM) scenario has proven to be very successful 
on large galactic and extragalactic scales, from cosmic str uctu re formation to the outer rotation 
curves of galaxies and the stability of galactic disks (Ostriker Sz Peebles 1973). The success of the 
CDM model however is also to some extent frustrating as any additional physical properties of the 
CDM particle remain hidden. Scientists are therefore searching for failures of CDM predictions 
that might lead to new insight into the nature and origin of the dark matter particle. 

Several small-scale problems of the CDM model have been discussed in the literature^ These 


inclu d e the distribution of satellite galaxies in large planar structures (e.g. iKroupa 


2005 


2013; 


Ibata et al 


Goerdt et al. 2014 ) or the mass-luminosity problem of satellite galaxies ( KrounaT 201flh . 


2011 


Qgiva fc Mori 2014 ; lOgiva et ah 2014 ; Qgiva fc Burkertl 2015 ). While simulations predict 


Flores V Primack 

1994; 

Burkert 1995: 

Strigari et al. 

2008 

Primack 2009; 

de Blok 2010; 

Bovlan-Kochlin et al. 


cuspy central density profiles of CDM halos with the density increasing steeply towa rds the center 


(IDubinski fe Carlberdll99ll : iNavarro et al.lll997l : iMoore et al 


often indicate a flat dark matter density core (e.g. 


1995 


de Blok et ah 


2008 


Gentile et al 


2009: 


1999; 


Dekel et al 


Flores &: Primack 


Oh et al 


1994 


Moore 


2003) , obs e rvations 


1994: Burkert 


2011). The detection of cored dark matter 


halos is not necessarily inconsistent with CDM. Various mechanisms have been identified that can 


in the galactic potential, induced by AGN feedback and galactic winds (e.g 

. Navarro et al 

1996; 

Ogiva & Mori 2011, 2014; Pontzen & Governato 

2012; 

Teyssier et al. 

2013; 

Amorisco et al. 

2014) 

or gaseous and stellar chimps spiraling to the center (e.g 

. El-Zant et al. 

2001; 

Via & Bovlan-Kolchin 

2004: Tonini, Lapi & Salucci 2006; Goerdt et al. 

2010; 

[none & Saitoh 

2011) 



The cusp-core problem is best documented in low-mass dwarf galaxies which are characterised 
by low baryon fractions and which therefore are ideal tracers of the underlying dark halo structure, 
unperturbed by the gravitational influence of the baryonic component. In addition, dwarf galaxies 
appear to host the highest density dark matter cores, making an an alysis of the halo structure 
easier. Ideally one would like to investi gate Ha/HI rotation curves (e.g. ICarignan fc Freemanlll985l : 
de Blok et al.ll200ll : ICorbelli et al.l I2014J ) which are a clear tracer of the gravitational potential as 
function of radius. One of the problems however is that for low-mass galaxies the stellar- and 
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gas velocity dispersions begin to exceed their rotational velocity. This is especially true for one 
of the smallest galaxies, known as dwarf spheroidals (dSph), the target of this paper. Recently, 


Kormendv fc Freemanl (KF14; 2014) investi gated the_dark halo scaling laws in late-type galaxies, 
including dSphs (see also Burkert 1997; Salucci et ah 2012). Interestingly, the stellar velocity 
dispersion of the dSphs is of order 8-10 km/s, very similar to the universal t urbulent velocity 


of the diffuse gas component in most low-redshift star forming disk galaxies (Dib et al. 2006). 


However, in contrast to more massive galaxies, the gravitational field in dSphs is small, generating 
rotation curves of order v ro t ~ 10 km/s (see section 4). Even if the star forming gas would have 
had enough angular momentum to settle initially into a thin disk configuration, stellar feedback 


processes could easily destroy the disk,_e 


spheroidal stellar system (jNavarro et al 


ecting gas in a wind and leading to a dispersion-d ominated 


1996; 


Mailer &; Dekell 120021 : iTevssier et al.l 120131 ) . 


KF14 consider a spherically symmetric stellar system with constant <r*, embedded in a spheri¬ 
cally symmetric dark matter core with constant density po,d- Solving the hydrostatic equation (see 
equation 3) and assuming that the stellar system is not self-gravitating they find that the stellar 
density distribution should be a Gaussian 


p*(r ) = po ,* x exp 



with po,* = p*(r = 0) the central stellar density and 


CL* 


2,al \ 1/2 

2vrGp 0 ,dy 


( 1 ) 


(2) 


the scale length of the stellar system. The projected stellar surface density distribution £* in this 
case is also a Gaussian with the same scale length. Remarkably, KF14 show that many dSphs 
follow a Gaussian surface density distribution better than the typical exponential profile, seen for 
galactic disks. 

KF14 however also find that some dSphs cannot be fitted by a Gaussian. In addition, also 
dSphs with inner Gaussian slopes show deviations further out. Another problem is the fact that 
isothermal, isotropic dark matter cores in equilibrium cannot have precisely constant densities, 
requiring a more detailed investigation. Finally, equation 2 does not provide any information about 
the dark halo core radii, velocity dispersions and masses. KF14 shift the observed stellar scale 
length a* and velocity dispersion cr* along lines of constant po,* onto the core scaling relations of 
more massive galaxies, in order to infer the halo core properties. It is however not clear whether 
dSphs should follow the same core scaling relations as more massive galaxies. 

In this paper we therefore have a more detailed look at the coupled kinematics of stars and dark 
matter in dSphs, relaxing the assumption of a constant density dark matter core. Section 2 solves 
the hydrostatic equation of two isothermal particle systems, coupled by their joint gravitational 
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field. We derive formulae on how to determine the dark matter core density and how to shift the 
observed stellar velocity dispersion c* and central stellar scale length a* in order to infer the dark 
halo velocity dispersion ad and halo core radius r C} d- Section 3 then focusses on deviations from 
Gaussian profiles and the origin of King profiles in isothermal dSphs. This section demonstrates that 
the stellar King concentration parameter c* is tightly related to the dark halo velocity dispersion 
(Jd and by this also makes it possible to determine the halo core radius and core mass. These 
analytical results are then applied to the Milky Way’s system of dSphs in sect ion 4_to investigate 
their halo core properties. Section 5 discusses the conjecture of Strigari et ajJ (2008!) that dSphs 
have a universal mass, of order 10 7 M 0 , within a scale radius of 300 pc. Section 6 summarises the 
results and concludes. 


2. The structure of double-isothermal (DIS) particle systems 

We assume that the cores of dark matter halos in dSphs are isothermal and isotropic with a 
constant velocity dispersion ctq^. This assumption certainly has to break down at some point as 
otherwise the dark halo mass would diverge as M ( }{r) ~ r. For investigations of the outer halo 
regions, the Burkert profile ([Burkertl 19951) might therefore provide a better approximation. It 
combines an isothermal-like inn er core with the characteristic W 3 outer density decline, seen in 
most CDM simulations ({Navarro et al.l 19971 ). 


Observations also indicate that the stellar body of dSp hs is characterised by an almost con¬ 
stant velocity disp ersion <7* , well beyond the half-l ight radius (jWalker et al.ll2009l : lEvans et al.l l2009; 
Salucci et al .11 2 012 . KF14). iMaiewski et al.l (|2013l ) recently reported a dro p in a* in the heart of th e 
Sagittarius dwarf spheroidal galaxy. A similar feature is seen in Sculptor (Bredd els fe I lelm i 20141), 
indicating that some dSphs might be more complex 2-component galaxies (lAmorisco fc Evans 
2011). The dense cores are limited to the very center that is small compared the half-light ra¬ 
dius. A central decline in c* might in fact be a characteristic property of dSphs in general. Many 
of them show a strong increase in stellar surface density in the very center, requiring a change 
in their kinematics, most likely a declining velocity dispersion, in order to be in hydrostatic equi¬ 
librium. The origin of cold nuclei in dSphs is certainly an interesting and yet unsolved problem. 
Here, however, we are interested in the global structure of dSphs which to good approximation 
is observed to be isothermal and will neglect their cold hearts. In addition, in order to keep the 
number of free parameters to a minimum, we assume that both, the stellar component^ and the 
dark matter cores are isotropic with negligible anisotropy effects (ICiottil 11993 : Evans et al.l 12009 : 
Salucci et al. 2012). 


Adopting a central density pd(r = 0) = po,d> the radial density distribution pd(r) of an isother¬ 
mal halo in hydrostatic equilibrium is determined by 


2 dlnp d 

° d dr 


GM d (r) 


( 3 ) 
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Here a d is the dark matter velocity dispersion and M d {r) is the cumulative dark matter mass inside 
radius r. The assumptio n of virial equilibrium might not always be va lid (IKroupalll997l ). especially 
for tidal dwarf galaxies (|Wetzstein et al.ll2007l: iPlo eckinge r et al.ll2015ll and s trongly tidally af fected 
dSphs like the Sagittarius dwarf spheroidal (jlbata et al.lll994J : iKroupall 19971 : lYang et al .1120141 ). De¬ 
viations from hydrostatic equilibrium might therefore help to identify those (mostly outer) regions 
of dwarf satellites, that are tidally interacting with the host galaxy. The thick grey line in the upper 
left and right panels of figure 1 shows the density distribution and logarithmic slope din pd/dlnr , 
respectively, of a non-singular isothermal sphere with finite po,d- Inside the core radius that for an 
isothermal sphere is defined as 





(4) 


the dark matter density distribution is roughly flat. Note however that Pd(r) is not exactly constant. 
This is due to the fact that, in equilibrium, the pressure gradient a^dln(pd)/dr has to balance the 
gravitational force GMd(r)/r 2 . One might consider this a negligible effect for r < r Ct d- It is however 
precisely this small density gradient that determines the density distribution of the embedded stellar 
component. 


Let us therefore now include a non self-gravitating, isothermal stellar component, as observed 
for dSphs. The stars then represent direct tracers of the underlying dark matter potential. It 
is unlikely that dSph s started that way. Obviously, the gas clouds from which the stars formed, 
were self-gravitating (jBurkert, fc Hartmannl 120131 : iNipoti fc Binnevl 120141 ) . A low star formation 
efficiency, combined with strong galactic winds and ram pressure stripping might however have 
removed most of the baryons before th ey cou ld condense into stars. This conclusion is also supported 
by the low nretallicites of dSphs (Dejkel & Silk 1986). If the galactic outflow was violent enough, 
in addition to leaving behind a non self-gravitat ing stellar system, i t could also have res haped an 
initially cuspy dark halo, generating a core (e.g. iNavarro et al.lll996l : iTevssier et al .1120131 ) . 


As both, the stars and the dark matter particles move inside the same joint gravitational 
potential their density distributions are coupled: 


2 d In p* GM d (r) 2 dln Pd 

= —— = <,d ~dr’ (5) 

Here p*(r) and <r* correspond to the density distribution and velocity dispersion of the stellar 
component. Equation 5 leads to 


P*( r ) =4x p K d (r ) 


( 6 ) 


with 
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(7) 


A is a constant of integration that determines the total stellar mass. Equation 6 shows that it is 
indeed the dark matter density gradient that determines p*(r ): 


dlnp* d In p d 

~T\ - = -' ( 8 

d In r d In r 

The three dashed lines in the upper panels of figure 1 show p*(r) for embedded stellar systems of 
our DIS model with velocity dispersions cr*/(Jd of 0.7, 0.5 and 0.25, corresponding to values k of 2, 
4 and 16, respectively. For k > 1.5 the stellar density distribution decreases faster than r~ 3 in the 
outer region. In the academic limit that the DIS model holds for all r, the stellar system would 
have a finite mass, despite the fact that it is isothermal at all radii. 

The middle, left panel of figure 1 shows the logarithm of the stellar density profile p*(r) and 
surface density distribution E*(r) (dashed and solid lines, respectively) as function of r 2 for k = 2,4, 
and 16. A Gaussian profile would be represented by a straight line and indeed fits the inner profiles 
in general quite well. The larger k. the similar are the density and surface density distributions 
and the more do the profiles resemble a Gaussian. 

Despite the fact that the dark matter core density is not precisely constant and by that 
equation 2 is not exactly valid we can still formally derive an approximate Gaussian scale length 
a* = — (din E*/cir 2 ) -1 / 2 by a least-squares linear fit of InE* versus r 2 within the innermost regions 
of the stellar component that we define as region, where E* decreases by a factor e with respect 
to the central value. The middle right panel of figure 1 shows that a*/r Cj ^ depends strongly on k. 
Solving for p 0 ,d i n equation 2 and inserting it into equation 4 we find 



It is not clear whether this should work, given the fact that equation 2 was derived for a con¬ 
stant density core while we have argued that it is actually the dark matter density gradient that 
determines the stellar density distribution (equation 8). However the solid black line in figure 1 
shows that equation 9 indeed provides an excellent fit to the actual data, derived from a numerical 
integration of the DIS model (red points). 

In the lower left panel of figure 1 we test the validity of equation (2) as an estimation for the 
underlying dark halo density po,d, given a* and cr*. If we write 


Po,d — Vp 


3u 2 

27rGa 2 


( 10 ) 
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log r/r cd log r/r c-d 




( r / r c,d ) 2 


log (a,/r cd ) 



K 


Fig. 1.- The thick lines in the upper left and right panels show the normalized density distribution 
and logarithmic density slope of an isothermal dark halo as function of radius. The dashed lines 
depict three embedded stellar systems with a*/<Jd = 0.7, 0.5 and 0.25, respectively. The thick line 
and the solid lines in the middle left panel show the surface density distribution of the dark halo and 
its three embedded stellar systems, respectively. Dashed lines in the middle left panel correspond 
to the normalized density distribution which is very similar to the surface density profile, especially 
for cold stellar systems. The red points in the middle, right panel show the correlation between 
c7*/f7rf of dark matter confined stellar systems as function of a*/r Ci d, as determined by integrating 
the equations 3 and 6. The solid line shows equation 9. In the lower left panel the red dots show 
the correlation between the dark matter density parameter rj p (equation 10) and k, derived by an 
integration of the DIS model. The solid line represents the empirical fit formula, equation 11. 











an analyses of the DIS model shows that rj p depends only on k and in a way as shown by the solid 
black line. For s>la very good approximation is (red points) 

rip « 1.01(1 + 0.5exp[— (k, — l) 0 ' 6 ]). (11) 

As expected, for very cold stellar systems with k> 1 the stellar system traces the innermost dark 
halo core with an almost constant density distribution. Here r/ p ~ 1 and equation 2 provides a good 
estimate of po,d- For kinematically hotter stellar systems, however, equation 2 is not valid anymore 
and the correction factor r] p has to be taken into account. 


3. King profiles and apparent extra-tidal components 


Up to now we focussed on the innermost regions of dSphs that are sensitive to the central dark 
matter density. In order to gain information about the dark halo core radii and core masses we need 
to look at stellar traces further out. According to equation 6 very cold dSphs with k> 1 populate 
regions that are deeply embedded in the dark halo core and that are therefore not good probes 
to explore the larger environment. One such example is Carina, shown in the upper left panel of 
figure 2. Carina can be fitted well by a Gaussian over most of the stellar body with a small change 
in slope dln£*/dr 2 in the outermost region. In contrast, the stellar distribution of the much more 
extended Sculptor dSph (upper right panel of figure 2) deviates strongly from a Gaussian. The red 
solid lines in both figures show the surface density profiles of DIS systems following equations 3 to 
7 with k, = 3.3 and k = 1.5 for Carina and Sculptor, respectively. With respect to the dark matter 
component, the stellar system in Sculptor (a* = 0.81fj,i) is kinematically hotter than in Carina 
(a* = 0.55arf), leading to a more extended structure. 


It turns out that the equations 3 to 7 lead to surface density distributions that fit all dSphs very 
well, even those with more complex, extended stellar components like_Sculptor. This results_from 
the fact that dSphs are in general observed to follow King profiles (King 19661; A morisco &; Evans 


2011) with different concentration parameters c*. King profiles are a 1-parameter family, charac¬ 


terised by the concentration parameter c which is equal to the logarithm of the ratio of the outer 
cutoff radius to the core radius of a particle system. The same is true for the projected surface 
density distributions of stellar systems in our DIS model. As an example, the solid lines in the left 
panel of figure 3 show three different stellar systems of our DIS model with n values of 2.8, 1.4 and 
1, respectively. The points show the corresponding best fitting King profiles which have concentra¬ 
tions of 0.6, 1.25 and 3.6, respectively. Note that this excellent fit hides a fundamental difference 
between DIS models and King models. Our stellar systems have constant velocity dispersions. 
They are part of a 2-component system, with the stars being embedded in a surrounding dark halo 
that has in general a different velocity dispersion than the stars. King models, instead, are one- 
component, self-gravitating particle systems that are sometimes also called truncated isothermal 
spheres. They are characterised by a special velocity distribution function that has been designed 
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r 2 [arcmin 2 ] r 2 [arcmin 2 ] 


Fig. 2.— They blue points in each panel show the stellar surface density measurements of 
Irwin & Hatzidimitriou (199a) of the Milky Way dSphs Carina and Sculptor. DIS models (red 


lines), normalized to the observed central surface density and scale length, with k values as given 
by table 1 provide an excellent fit to the data. 
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to fit stellar systems like globular clusters with sharp outer edges, generated as a result of tidal 
stripping. In order for such a sharp outer edge to exist, the velocity dispersion of the stars in the 
King model has to decrease with radius with <r* —y 0 at the outer edge. Otherwise stars would 
be able to move beyond it. It is therefore surprising and at first not necessarily expected that the 
one-component King profiles with a completely different kinematics provide such a good fit to the 
stellar structure of our two-component DIS systems, over more than 4 orders of magnitude in In X*. 

A characteristic property of King models is that the surface density structure changes strongly 
for concentrations 1.2 < c* < 2 from a core with a steeply decreasing outer edge to a more 
extended structure. The DIS models follow this trend nicely with kinematically hotter stellar 
systems, characterised by smaller values of ft, corresponding to King models with larger c*. In 
the transition regime, however, the best fitting King profiles are somewhat steeper than the stellar 
systems for lnX*/Xo,* < —4 (see the c* = 1.25 profile in the left panel of figure 3). Interpreting 
such an extended population of stars as extra-tidal, in this case, would be misleading. These stars 
are still deeply embedded and strongly bound to their dark halo. On the other hand, an extra-tidal 
component detected at that level in lnX*/Xo,* in systems with concentrations c* < 1 or c* > 2 
cannot be explained within the framework of our model and therefore might indeed represent a 
separate hot halo or even an extra-tidal component. 

The right panel of figure 3 shows that there exists a tight correlation between the k value of the 
DIS model and the King concentration parameter c* (solid black line). The kinematically hotter 
the stellar system, i.e. the smaller ft, the more extended the stellar system and the larger c*. The 
red points show the empirical relation 


log ft = 1.25 exp(—1.72c*) (12) 

which is an excellent fit to the data. Given a*, er* and c* one can now use the equations 9 to 12 
and calculate the central density po.d; velocity dispersion and core radius r c ^ of the dark halo: 


log 


log ^ 
V P c 

<?d 

km/s 

Po,d 

M 0 /pc 3 


log 


0.088 + log + 0.625exp(—1.72c*) 


= 0.625 + log 


= 2.04+ 2 log 


O * 


km/s 
o * 

km/s 


+ 0.625 exp(—1.72c*) 


(13) 


- 2 log 


(s) 


4. The dark halo core properties of local Milky Way dSphs 


As an application we investig ate the dark halo c ore properties of the 8 classical Milky Way 
dSphs (table 1), observed by Irwin & Hatzidimitriou (1995, IH95, see also Kormendy & Freeman 
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In r/r cd c * 

Fig. 3.— The left panel compares the normalized surface density distribution of three characteristic 
King profiles (red points) with different concentrations c* = 0.6 (innermost contour), 1.25 (middle) 
and 3.6 (outer contour) with the stellar surface density distribution of our DIS model (solid lines). 
Note the apparent ’’extra-tidal” stellar component for the profile with c* = 1.25. The solid line in 
the right panel shows the correlation between the k values of the DIS model and the concentration 
of the best fitting King profile. The red points show the empirical fit formula, equation 12. 


Table 1: Physical properties of the stellar and dark halo component of the 8 classical Milky Way 
dSphs (KF14, IH95, McConnachie 2012). 


stellar component dark matter component 



a* [pc] 

cr* [krn/s] 

c* 

r c ,d [pc] 

a d [krn/s] 

Po,d [Mq/pc 3 ] 

M Ctd [10 7 M 0 ] 

M 300 [ 10 7 M 0 ] 

Carina 

202 

6.6 

0.51 

450 

12.0 

0.13 

2.6 

0.9 

Draco 

176 

9.1 

0.50 

397 

16.7 

0.33 

4.4 

2.1 

Leo I 

221 

9.2 

0.58 

460 

15.6 

0.22 

4.6 

1.6 

Leo II 

174 

6.6 

0.48 

400 

12.4 

0.18 

2.4 

1.2 

UMi 

211 

9.5 

0.51 

470 

17.3 

0.25 

5.6 

1.9 

Fornax 

705 

11.7 

0.72 

1310 

17.7 

0.035 

17.3 

0.4 

Sextans 

400 

7.9 

0.98 

641 

10.3 

0.053 

3.0 

0.5 

Sculptor 

189 

9.2 

1.12 

286 

11.4 

0.33 

1.7 

1.5 
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2014). IH95 determined their stellar surface density distribution with high enough resolution in 
order to derive King concentration parameters and determine the central scale lengths. IH95 also 
provide stellar velocity dispersions. The central Gaussian scale lengths a ma j were determined from 
major axis surface brightness profiles, shown in figure 2 of IH95. These values are very close to 
the King core radii, summarized in Table 4 of IH95. Following IH95, a* was then derived as the 
geometrical mean along the major and minor axis with a* = a ma j x r c , g jH/r c jH where r c%g jH and 
t c jh is the stellar system’s geometric mean and major axis core radius, respectively, as determined 
by IH95. 


Using the set of equations (13) we now can calculate the halo core parameters r Ct d, &d and po,d- 
The results are summarized in table 1. The upper left panel of figure 4 shows <7* versus a* (grey 
triangles) and versus r c ^ (red points with error bars) for the 8 dSphs. Since IH95, u pda ted 
stella r velocity dispersion measurements have been published (e.g. lWalker et al.ll2009l : lMcConnachie 


2012). Here we ado pt the values of er* as given in the regularly updated McConnochie data base 


(McConnachie 2012). Dark halo cores are hotter than their stellar systems with velocity dispersions 
in between 10-18 krn/s and, on average, <7d ~ 1.6cr*, correspondin g to k ~ 2.56. The halo core radii 


lie in the range of 290 pc to 1.3 kpc. On average, r C} d ~ 1.9a*. D onato et al. (120041 ) analysed a 


sample of high-resolution rotation curves of 25 disk galaxies and determined independently the disk 
scale lengths rdisk and dark matter core radii. They found that both radii are strongly correlated 
with rdisk ~ 2.4 r C} d- It is remarkable that the ratio between the stellar scale length and the dark 
halo scale length is the same (of order 2) in very different galactic systems over a large range of 
mass. The origin is still unclear and might provide further insight into the mechanisms that lead 
to dark matter cores. 


The upper right panel of figure 4 shows the dark halo central surface densities po,d as function 
of r C) rf. Typical values are 0.2 M 0 pc~ 3 with a range of 0.03 - 0.3 M 0 pc -3 . It has been argued 


(e.g. 

Athanassoula et al. 

1987; 

Burkert 

1995; 

Salucci & Burkert 

2000 ; 

Kormendv & Freeman 

2004; 

de Blok et al. 

2008; 

Gentile et al. 2009; 

Donato et al. 2009 

; Kormendv & Freeman 2014 

see however Saburova & d 

2014 

; Cardone & del Popolo 

2012). KF14 find that this scaling relation holds over more than 18 


magnitudes in Mr. The black triangles in the upper right panel and in both lower panels of figure 4 
show the core properties of galaxies, compiled from the literature by KF14 (see the list of references 
for the original data in table 1 of KF14). The core surface densities of all galaxies lie in a narrow 
range of 30 M 0 /pc 2 < (po,d x r c ,d) < 160 M 0 /pc 2 (dashed lines). The 8 Milky Way dSphs fall 
precisely into this regime, despite the fact that their r c d are on average a factor 6 smaller with p$d 
being a factor of 6 larger. The lower left panel of figure 4 shows the dSph core masses, which for 
non-singular isothermal spheres are 


M c .d = 2.17 x po,d r 3 c4 = 162.75 


/ {po,d x r Cjd ) \ 
V 75M 0 pc- 2 ) 



(14) 


The core masses cover a range of one order of magnitude with masses in between (see table 1) 
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Fig. 4.— The grey triangles in the upper left panel show the stellar velocity dispersion <r* versus 
the central stellar scale length a* of the 8 classical Milky Way dSphs. Red points with errorbars 
depict the corresponding dark halo velocity dispersion versus the halo core radius r C: d- The 
red points in the upper right and lower left and right panels show the correlation of r Cy d with the 
dark halo central densities po,d> core masses M c and velocity dispersions respectively. Black 
triangles depict the core properties of more massive galaxies. Typical error bars for these data 
points are shown in the upper right or left corners of each plot. The thick grey line in all four panels 
corresponds to the dark halo core scaling relation (/9o,<Rc,d) = 75 M 0 pc -2 that fits all galaxies very 
well. The two dashed lines show the observed spread with the upper and lower limit corresponding 
to core surface densities of 30 and 160 M 0 pc~ 2 , respectively. 
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1.7 x 10' < M c( j < 1.7 x 10 8 M 0 . The expected correlation between M cr i and r c ,d (equation 14) is 
drawn for (pod x r c,d) = 30, 75 and 160 M 0 pc” 2 , together with the core masses of more massive 
galaxies. The dSphs follow the same core mass scaling relations as massive galaxies with the same 
spread. Finally, the lower right panel shows again ad versus r Ct d- Now we compare the dSphs with 
the more massive gala xies. Both follow the same universal scaling relation cr ( 2 x r”^ = 0 . 45+0 27 
(km/s ) 2 pc ” 1 ([see also de Vega et al. 120141. again with precisely the same spread. 


5. The origin of a common mass and length scale for dark matter cores 


S tr igari et ah (12008 ) proposed that all dSphs of the Milky Way have the same total dark matter 
mass log(M 300t d/M & ) = 7.0 Iq 8 , contained within a radius of 300 pc. The origin of a universal 
and constant mass might at first appear surprising, given the fact that the core masses are a strong 
function of core size (equation 14). Note however that M c j is measured within r c j whereas M 3 oo,d 
is the mass within a fixed radius 300 pc that can be smaller or larger than r C ( i- The question still 
arises why there should exist such a universal radius r u d, inside which halo cores have the same 


mass A/,, d and what determines this radius. In addition, adding Andromeda dSphs, Collins et al 


(120141 ) find outliers that are not consistent with the Strigari et al. mass which indicates that the 
situation can be more complex. 


Qgiva et al.1 (|2014l ) discussed a possible connection between t he existence o f a un iversal mass 


scale and the universal core surface density of dSphs. Following Ogiva et a l. (120141 ). let us now 
explore this question within the context of the DIS model. We start with a population of dSphs 
that has a common core surface density (po,d x r c ,d). Given r c ,d, one can determine po,d and with 
equation (4) ad ■ An integration of equation (3) then gives the complete density profile and by this 
M v d for any given value of r U: d- The horizontal lines in figure 5 show Af 3 oo,d as function of r C: d for 
the typical surface densities (po,d x fc,d) = 75+ 4 q M 0 pc ” 2 of the Milky Way dSphs. Interestingly, 
M 300 ,d is not continuously increasing with r Ct d but instead has a maximum at the point where it 
crosses the correlation between halo core mass and core radius (vertical lines). This point also 
corresponds to the adopted characteristic scale length r Ut d = 300 pc. The blue and red points 
in figure 5 show the dSph’s Af 30 o,d and M Cj d, respectively. We added 10 dSphs from table 1 of 
Strigari et al. (2008|) with given M 3 oo,d and King radii Thing ,* (blue triangles). The halo core radii 
are not known for this sample. For our joint sample we find on average r c ,d ~ 2.3 Thing,* which 
was used in figure 5 in order to estimate r c ,d for the dSphs with unknown core radii. The red 
points follow the expected correlation between core mass and core radius (vertical solid and dotted 
lines). The A/300,<i masses however show no such correlation with r c ,d but instead are roughly 
constant. Table 1 summarizes the A/300, d values of our sample of Milky W ay dS phs. With on 
average A/ 300 ,d = 1.3 ± 0.6 x 10” M 0 they are in excellent agreement with Strigari et al. (12008 ). 


The horizontal lines in figure 5 however show also that A/ 300 ,d is n °t precisely constant but 
should depend on the halo core radius. This is true for any adopted radius r U) d■ Within the 
framework of the DIS model there exists no universal radius r U} d inside which M u g, is constant for 
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a population of dSphs with given core surface density. It is however interesting that M ud reaches 
a maximum for halos with core radii r Cj d = f' u d . This is true for any adopted length scale r u d- As 
M u : j is very insensitive to r c ^ in the flat vicinity of this maximum, all dSphs with core radii in this 
regime would show very similar M Uy d values, which is exactly what we observe for the Milky Way’s 
dSphs. The best choice of r u ^ therefore is the average logarithmic core radius of a given sample 
of dSphs with universal core surface densities. M u> d is then the mass within a dark halo core with 
r c ,d = r u ,d (equation 14): 


log r Utd 


1 

N 


N 

^2 log r cAi 
1=1 


M U! d — 2.17(po, d x T c ^d)f u ^d 


(15) 


with N the number of dSphs and r c _d,i the core radius of galaxy i. The origin of a maximum for 
v Cj d = r U: d can be easily understood. For r c ,d > r u ,d the core radius is larger than the region sampled 
by r Ut d and the density is roughly constant with a value po.d ~ l/ r c,d due to the assumption of 
a constant core surface density. The enclosed mass with a given fixed radius r Ut d is then M Ut d ~ 
Po,d ~ l/ r c,d leading to M U)d decreasing with increasing r C)d . For r Ct d < r u ^ the region is larger than 
the core and extends out to radii where the dark matter density distribution begins to decrease. 
To first order we can approximate M U)d now as the mass of a constant density core pd(r) = Po,d for 
r < r C) d plus the mass of an envelope with a power-law density distribution pd = po,d x ( fc,d/f ) 2 - 
As po,d ~ l/ r c,d we get M U) d ~ T^ d {r U) d/r Ct d — 2/3) which is a continuously increasing function with 
increasing r Ci ^. 


Applying equation 15 to our sample of dSphs we find r Uj d = 400 ± 100 pc and with (po,d r c,d) 
= 75 Mq pc~ 2 get M u d = 2.4 ± 1.4 x 10' M 0 , in good agreement with Strigari et al. ( 2008 ). If the 
population has a large spread in r c ^ we also expect to find outliers, populating the wings of the 
distribution further away from the maximum with smaller masses M„ d- This might explain the 
detection of outliers in Andromeda’s system of dSphs ( Collins et al. 20141). 


6. Summary and conclusions 


Motivated by the conjecture that dSphs are isothermal stellar systems (e.g. Evans et ajJ 2003 ). 
embedded in isothermal dark matter cores we investigated the structure of two particle systems 
with constant but different velocity dispersions, in virial equilibrium within their joint gravitational 
potential. Note that here the main objective was not to demonstrate that dark m atter cores ar e 
isothermal and isotropic and it certainly has to break down outside of some radius (Burkert 19951 ). 
We worked with this assumption because it is the most simple model of a halo core with the least 
number of free parameters. Of course, the fact that we find solutions that fit the observations 
very well is promising. But it is not a proof due to the fact that a fine tuned radial distribution 
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log r Cid [kpc] 


Fig. 5.— Red and blue points show the dSph’s core masses M c ^ and masses M 30 o,d within a hxed 
radiu s r„ ,i = 300 ptxsummarized in table 1. Blue triangles show additional dSphs taken from table 
1 of Strigari et ah (2008). The red points follow the core scaling relations po t d x r C: d = 75^40 M P 


■0 


pc~ 2 , observed for massive galaxies (solid and dotted lines). The second set of solid and dashed 
lines with a maximum at r C: d = 300 pc shows the predicted relationship between M 30 o,d and r Ct d 
for halos with the same core surface densities. All lines are labeled according to their core surface 
density in units of [M 0 pc~ 2 ]. 
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of anisotropy, coupled with a properly chosen gradient in velocity dispersion could always lead to 
similar cored density profiles. 


We demonstrate that the surface density distributions of the non self-gravitating stellar com¬ 
ponent of a DIS galaxy can show a rich variety of profiles. They can formally be fitted by King 
profiles despite the fact that the DIS model consists of two isothermal components in contrast to 
the one-component, non-isothermal King model. We find that the stellar systems in the DIS models 
have steeply decreasing outer edges, especially for high values of k, not because they are tidally 
limited but because they are deeply embedded within the inner core regions of their dark halo. 
This is in contrast to the real one-component King model which is characterised by a tidal radius, 
the Roche radius, where the gravitational potential of the host galaxy begins to dominate and 
where stars are unbound to the satellite. In addition, DIS systems have projected stellar velocity 
dispersion profiles that remain constant all the way to their outermost radius. One-component 
King systems, on the other hand, show outer velocity dispersion profiles that decrease approaching 
<7* = 0 at the tidal radius. Measurements of the stellar velocity dispersion of dSphs close to their 
cutoff radius therefore could help to distinguish tidally truncated one-component systems without 
confining dark halos ( Yang et ah 2014]) from those where the maximum radius is determined by 
a strong dark matter confinement. This is also important, as the outer radii of satellite galaxies, 
interpreted as tidal radii, have been used in order to gain information about the satel li tes’ orbital 
parameters or the dark halo mass distribution of the host galaxy (e.g. Pa setto et al. 2011). Our 


analyses instead shows that the maximum radii of the satellites measure the core radii of their own 
dark halos and therefore do not contain information about the tidal radius and the strength of the 
host galaxie’s tidal field. 


That ry* traces the dark halo’s core radii follows from the fact that the dSphs on average have 
k ss 2.6. ry* is close to the point, where the logarithmic density gradient dinp*/dinr begins to 
decrease faster then -3 which, according to equation 8, then corresponds to a dark halo density 
gradient of dlnp/dlnr = —3/k = —1.2. The upper right panel of figure 1 shows that this slope is 
close to the core radius r c of the dark matter halo. It is not clear yet, whether this is a coincidence 
or whether the processes that generated the dark halo c ores and their non-self gravitating stellar 
tracer component naturally lead to such a configuration (jDekel fc Woo 2003). 


For DIS systems, the central Gaussian scale length a*, the velocity dispersion er* and the 
concentration c of the stellar component completely specify the dark halo core parameters />o,d; a d 
and r c d- We determined these parameters for 8 dSphs of the Milky Way and find that their dark 
halos have the same core surface densities /?o,d x r c,d = 75^45 M 0 pc -2 as more massive galaxies, 
with exactly the same spread. This is very puzzling as dSphs have a different structure and live 
in very different environments. At the moment it is not clear whether this result is true also 
for M31’s dSphs and whether it can be extended to the ultra-faint satellite population which due 
to their smaller radii should have even higher dark halo core densities. An analyses similar to 
what was presented in this paper would require deeper observations of their stellar surface density 
distributions that are accurate enough in order to make King profile fits. 
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itational interaction with the barvonic component (e.g. Navarro et al. 1996; El-Zant et al. 

2001; 

Goerdt et al. 2010; 1 

none & Saitoh 2011; Ogiva & Mor: 

2011; 

Pontzen & Governato 

20121: Governato et al. 

2012 

; Tevssier et al 

2013 

; Gritschneder & Lin 2013; 

Ogiva & Mori 

20 lj; 

Ogiva et al. 20L 

) to a 

non-standard primordial 

power spectrum 

(Zentner & Bullock 2002 

; Polisensky & Ricotti 

2014), 

warm dark matter (e.g. 

Lovell et al. 

2014 

), other intrinsic properties of dark matter like self- 

interaction and self-annihilation (e.g. 

Spergel & Steinhardt 

2000; I 

lurkert 

2000^ 

Loeb & Weiner 


2011 


Elbert et al. 2014 1 or modifications of Newtonian dynamics (e.g. MilgroirJ 


1983 


Kroupa 


2012). Whatever the mechanisms, considerable fine tune is required in order to generate a uni¬ 


versal core scaling relation over more than 18 orders of magnitudes in blue magnitude Mg with 
exactly the same spread. 

Adopting a constant core surface density, M c ^ depends strongly on r c . Focussing however on 
a fixed radius r Ut d, the enclosed mass M u ^ shows a different dependence on halo core radius, with a 
maximum at r c ^ = r u ^- All halos with core radii in the vicinity of this maximum should therefore 
show similar values of M u ^ which could explain the observations of Strigari et al. (.2008). The best 
choice of r u ^ is therefore dependent on the dSph’s distribution of r c ^. There does however not 
exist a universal mass scale M Ut d that is independent of r c ^. Smaller values of M u ^ are expected 
for outliers with core radii that are very different from r u Turning this argument around, if such 
a universal mass would exist, independent of r Cj d, it would be a clear signature that dark halo cores 
are not isothermal. 


As the core radii of dSphs are small, their core densities have to be high in order for the core 
surface density to remain constant. This should shield dSphs efficiently against the tidal forces of 
their host galaxies. Adopting a constant rotation curve v ro t , the Milky Way’s mean density within 
a given radius r is 


( Pmw) = 


3v 2 


rot 


4-7T Gr 2 


= 2.7 ( 


^ y/kpcy c _ 3 

y 220km/s J \r J 


(16) 


The stellar system in dSphs would be tidally affected if ( Pmw) > Po,d ~ 0.2 M 0 pc -3 which requires 
orbital pericenters of order a few kpc, which is very unlikely. dSphs therefore should be strongly 
shielded from any tidal affects by their deep dark matter potential wells and should survive as 
satellites of the Milk y Way for a long time to come. However extra-tidal debris has been reported 
in some dSphs (IH95 Walcher 2003 : Maiewski et al. 2005 : Battaglia et al. 2012 1. We discussed in 
section 3 that the DIS model in a certain concentration regime indeed leads to profiles that are 
somewhat more extended than the best fitting King profiles. This could be mis-interpreted as a 
tidal component. Strong evidence for tidal interactions would however represent a real c hallenge 
for the existence of a shielding dark halo, opening the door for alternative ideas (jMilgroml 19831; 
Yang et ah 2014). 
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